On the role of injection in kinetic approaches 
to nonhnear particle acceleration 
at non relativistic shocks 



S. Gabici*, P. Blasi^ and G. Vannoni* 

*Max-Plank-InstitutfilrKernphysik, Saupfercheckweg 1 , 691 17 Heidelberg, Germany 
' INAF/Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy 

Abstract. The dynamical reaction of the particles accelerated at a shock front by the first order 
Fermi process can be determined within kinetic models that account for both the hydrodynamics of 
the shocked fluid and the transport of the accelerated particles. These models predict the appearance 
of multiple solutions, all physically allowed. We discuss here the role of injection in selecting the 
real solution, in the framework of a simple phenomenological recipe, which is a variation of what is 
sometimes referred to as thermal leakage. In this context we show that multiple solutions basically 
disappear and when they are present they are limited to rather peculiar values of the parameters. 



Diffusive shock acceleration is thought to be responsible for acceleration of cosmic 
rays in several astrophysical environments. Despite the success of this theory, some is- 
sues are still subjects of much debate, for the theoretical and phenomenological impli- 
cations that they may have. One of the most important of these is the reaction of the 
accelerated particles on the shock: the violation of the test particle approximation oc- 
curs when the acceleration process becomes sufficiently efficient that the pressure of 
the accelerated particles is comparable with the incoming gas kinetic pressure. Both the 
spectrum of the particles and the structure of the shock are changed by this phenomenon, 
which is therefore intrinsically nonlinear (Ellison, these proceedings). Nonlinear effects 
in shock acceleration of thermal particles result in the appearance of multiple solutions 
in certain regions of the parameter space. This phenomenon is very general and was 
found in both the two-fluid |4] and kinetic models 10, SS- Here we investigate the 
phenomenon of multiple solutions and show that the appearance of these solutions is 
dramatically reduced if a self consistent model for injection is adopted. 



A SEMI-ANALYTICAL APPROACH TO THE PROBLEM 



Following the approach presented in , we solve the steady-state transport equa- 

tion for the cosmic ray distribution function f{x, p) at a plane shock wave: 
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coupled with the continuity and Euler equations describing the dynamics of the flow: 

pQUQ = pu ; pQul-\-Pg,() = pu^ +Pg+PcR. (2) 

All the symbols have their usual meanings, and the subscript refers to quantities 
measured upstream infinity. Pg and Pcr represent the contributions to the total pressure 
of thermal gas and cosmic rays. We introduce the quantity Up defined as: 

Up = ui-^ I dx^f{x,p), (3) 
Jo j-oo ax 

whose physical meaning is instrumental to understand the nonlinear reaction of particles. 
The function Up is the average fluid velocity experienced by particles with momentum 
p while diffusing upstream away from the shock surface. In other words, the effect of 
the average is that, instead of a constant speed ui upstream, a particle with momentum 
p experiences a spatially variable speed, due to the pressure of the accelerated particles. 
Since the diffusion coefficient is in general p-dependent, particles with different energies 
feel a different compression coefficient, higher at higher energies if, as expected, the 
diffusion coefficient is an increasing function of momentum. 

It can be shown that equations [T] and |2l can be reduced to an integral-differential 
equation for the quantity U (p) = Up/uQ, that when solved with the boundary condition 
U (pmax) = 1 ' provides us with the functions Up and foip) describing the flow profile 
and the particle distribution function at the shock |llll2.0]- 

In the problem described above there are several independent parameters. While the 
Mach number of the shock and the maximum momentum of the particles are fixed 
by the physical conditions in the environment, the injection momentum pi„j and the 



This corresponds to assuming that the fluid is not affected by the cosmic ray pressure at upstream infinity. 



acceleration efficiency 77 are free parameters. The procedure to be followed to determine 
the solution was defined in |2, 3]. In Fig. 2 we plot U (pmax) as a function of the total 
compression factor of the shock Rtot, for Tq = lO^K, pmax = lO^mc and pinj = lO^^mc 
in the left panel and pinj = lO^mc in the right panel (m here is the mass of protons). The 
parameter rj is lO^'^ in the left panel and 10^^ in the right panel. The different curves 
refer to different choices of the Mach number at upstream infinity. The physical solutions 
are those corresponding to the intersection points with the horizontal line U (pmax) = 1, 
so that multiple solutions occur for those values of the parameters for which there is more 
than one intersection with U (pmax) = 1 • These solutions are all physically acceptable, as 
far as the conservation of mass, momentum and energy are concerned. FigC] shows that 
multiple solutions are a very common phenomena, since they are found if very different 
values of the parameters are adopted. In the next section we show how the situation 
changes if a better description of injection is implemented in the model. 

A RECIPE FOR INJECTION 

The injection of particles into the cosmic ray population at a shock can be understood 
only considering the complex non-linear interactions between suprathermal particles, 
MHD waves and background thermal plasma. Due to its intrinsic complexity, the injec- 
tion process is often parametrized by means of an injection momentum pi„j, representing 
the minimum momentum of the particles that can be accelerated, and an efficiency rj , 
which fixes the fraction of the thermal particles that are injected in the accelerator. An- 
other possibility is to adopt the thermal leakage model to describe the injection LX^. In 
this model, the post shock gas is assumed to be thermalized at a temperature T2. Protons 
in the tail of the Maxwellian distribution can recross the shock and go back upstream if 
their velocity is high enough to allow them to avoid trapping by waves. Those protons 
are injected in the accelerator. Usually, in this model, the injection momentum is set to a 
few times the thermal momentum: pinj = ^ -JlmFTi with ^ tuned in order to fit numeri- 
cal 10] studies of diffusive shock acceleration. It is important to stress that the parameter 
^ fixes both the values of injection momentum and efficiency, which are no longer free 
parameters but are now connected in a physically motivated way. It is easy to implement 
such a recipe in the calculations L3J. 

RESULTS AND CONCLUSIONS 

The appearance of multiple solutions can be investigated in the whole parameter space, 
in order to define the regions where the phenomenon appears, when it does. In Fig. 8 
we highlight the regions where there are multiple solutions (dark regions) in a plane 
— log(pmax), for different values of the Mach number of the shock. The multiple 
solutions disappear in most of the parameter space, and when they appear they look 
as narrow strips in the parameter space, at the boundary between regions describing 
unmodified and modified shocks respectively [3]. This result may suggest that the 
narrow regions indicate the transition between two stable solutions, although this needs 
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FIGURE 2. Parameter space for multiple solutions. In the dark regions multiple solutions are still 
present. 



further confirmation tlirougli detailed analyses of the stability of the solutions. 
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